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PART  I. 


RECENT  ADVANCES  IN  MATHEMATICAL  MODELING  OF  THE  MECHANICS 
OF  BIMODULUS,  FIBER-REINFORCED  COMPOSITE  MATERIALS* 


1.  INTRODUCTION 

The  classical  elasticity  theory  hypothesis  [l]  that  the  material  obeys 
the  generalized  Hooke's  linear  stress-strain  relation  is  justified  for  appli- 
cation to  most  metallic  alloys  loaded  below  the  elastic  limit.  For  more  highly 
flexible  materials  including  polymers,  various  nonlinear  constitutive  relations 
such  as  neo-Hookean,  Mooney-Rivlin,  and  Ogden  have  been  proposed;  see  [2].  All 
of  these  relations  make  use  of  extension  ratio  (X)  rather  than  engineering 
strain  (e);  note  that  X - 1 + e.  Most  of  the  relations  have  a different 
elastic  modulus  at  a finite  value  of  compressive  strain  (cQ  < 0)  than  at  the 
same  absolute  value  of  tensile  strain  (|  e | > 0) . However,  in  agreement  with 
the  careful  measurements  of  various  investigators  (cf.  [3]),  these  finite-strain 
relations  have  no  discontinuity  in  slope  in  going  from  compressive  strain  (X  < 1) 
to  tensile  strain  (X  > 1) . 

In  contrast,  certain  materials  demonstrate  a distinct  change  in  modulus 
in  going  from  compression  to  tension*;  see  Fig.  1.  These  materials  apparently 
are  primarily  composite  materials,  as  listed  in  Table  1.  Jn  the  literature, 
this  class  of  materials  has  variously  been  called  bilinear,  blmodulus,  different- 
modulus,  and  multimodulus.  Here  the  term  bimodulus  is  believed  to  be  most  de- 
scriptive of  a material  having  different  linear  stress-strain  relations  in 
compression  than  in  tension. 

The  first  multidimensional  model  for  bimodulus  materials  was  proposed 
by  Ambartsumyan  [8]  for  Isotropic  material,  such  as  a composite  material  with 
spherical  particles.  It  was  later  extended  to  the  orthotropic  case  [9]. 

The  second  and  third  models  are  the  restricted  compliance  model  due  to 
Isabekyan  and  Khachatryan  [10]  and  the  first-invariant  model  of  Shapiro  [ll]. 

A fourth  model  is  the  weighted  compliance  theory  originated  by  Jones  [12J. 

The  fifth  model  is  the  fiber-governed  bimodulus  symmetric  compliance 
model  originated  by  Bert  [6]. 


In  the  next  section,  criteria  for  evaluating  bimodulus  material  models 
are  presented  and  in  subsequent  sections  the  criteria  are  applied  in  critically 
evaluating  the  various  models. 


★ 

This  part  of  the  report  is  a slight  expansion  of  a paper  of  the  same  title 
presented  at  the  15th  Annual  Meeting  of  the  Society  of  Engineering  Science, 
Gainesville,  Florida,  Dec.  4-6,  1978. 

Salnt-Venant  [4]  in  1864  made  perhaps  the  earliest  mention  of  material  with 
different  behavior  in  tension  and  compression. 


Tensile  Young's  Modulus 


Material 

Reinforcement 

Geometry 

Ref. 

Divided  by 

Comore.ssive  Young's  Modulus 

ATJ-S  graphite 

Granular 

5 

1.2 

ZTA  graphite 

Granular 

5 

0.8 

Glass-epoxy 

Fibrous 

5 

1.25 

Boron-epoxy 

Fibrous 

_ 

3 

0.8 

Graphite-epoxy 

Fibrous 

5 

1.4 

Carbon-carbon 

Fibrous 

5 

2.0  to  5.0 

Kevlar-rubber 

Fibrous 

6 

0.77  (transverse)  ^ 
to  305  (longitudinal) 

Polyester-rubber 

Fibrous 

0.75  (transverse)  A 

to  16.7  (longitudinal) 

Based  on  experimental  results  reported  by  Patel  et  al.  [7]. 


2.  CRITERIA  FOR  EVALUATING  MATERIAL  MODELS 


A linear  material  model  can  be  characterized  by  either  its  compliance 
matrix  or  alternatively  by  its  stiffness  matrix.  The  criteria  used  to  evaluate 
the  various  material  models  are  as  follows: 

1.  The  compliances  should  be  consistent  with  measured  values  for  the 
conditions  specified,  i.e.  they  should  depend  upon  the  nature  of  the  multlaxial 
stress  or  strain  state  as  appropriate. 

2.  The  shear  modulus  for  an  orthotropic  blmodulus  material  should  have 
different  values  for  shear  stresses  of  same  magnitude  but  opposite  sign  in  any 
coordinates  other  than  the  material-symmetry  directions. 

3.  It  would  be  preferable  to  have  the  strain  energy  be  positive 
definite  so  that  energy  is  conserved. 

Criterion  1 is  the  main  criterion  and  it  requires  that  the  model  be  able 
to  duplicate  measured  response  upon  changing  the  multlaxial  stress  or  strain 
state  as  appropriate. 

Criterion  2 is  clarified  by  Figures  2 and  3 due  to  Jones  [5].  In  Fig. 

2,  it  is  apparent  that  the  fibers  are  loaded  in  exactly  the  same  way  by  positive 
and  negative  shear  stresses,  while  in  Fig.  3,  it  is  clear  that  the  fibers  are 
loaded  in  tension  by  a positive  shear  stress  and  in  compression  by  a negative 
shear  stress. 

As  alluded  to  by  Voigt  [l3]  and  shown  by  Eubanks  and  Sternberg  [14]  and 
Lempriere  [15],  criterion  3 implies  that:  (a)  the  compliance  matrix  be  symmetric 
and  (b)  certain  limits  exist  on  the  compliances  so  that  the  compliance  matrix 
is  positive  definite.  Symmetry  of  the  compliance  matrix  is  necessary  in  order 
for  a material  to  be  mechanically  stable,  as  shown  by  Brun  [16].  Furthermore, 
compliance  symmetry  is  highly  desirable  in  that  most  structural  analysis  al- 
gorithms are  based  on  this  assumption,  i.e.,  they  have  no  provision  for  unsym- 
metrlc  compliance  or  stiffness  matrices. 

In  the  sections  to  follow  only  the  orthotropic  versions  of  the  various 
material  models  are  presented  and  these  are  limited  to  the  plane-stress  case  for 
brevity. 


3.  MODEL  I:  THE  AMBARTSUMYAN  MODEL 


In  this  theory  [9], 
as  follows 

G1  ■ bu  "l  + b12 


the  srains  are  expressed  in  terms  of  the  stresses 
°2 


C2  " b21  °1  + b22  °2 


(1) 


el2  " b61  °1  + b62  °2 


Here  , 02  are  principal  stresses;  e^,  e are  t*ie 
strains  associated  with  principal-stress  directions  1,  2 ; 
liance  coefficients  which  take  on  values  as  follows: 


if 


> 0 
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normal  and  shear 
bij  are  the  com- 
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Fig.  2 Positive  and  negative 
shear  stress  applied 
along  material -symmetry 
directions. 


Fig.  3 Positive  and  negative 
shear  stress  applied  at 
45  degrees  to  the  material- 
symmetry  directions. 
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In  a bimodulus  material,  b J ^ i b ^ by  definition,  while  in  a classical 
(-)  (+)  ^ 

material  b^  ■ b^  . It  is  noted  that  in  an  orthotropic  material  the 

principal-strain  directions  do  not  coincide  with  the  principal-stress  directions; 
thus,  the  presence  of  e^2  is  necessary. 

It  should  be  noted  that  the  lack  of  the  explicit  presence  of  a shear 

stress  does  not  mean  that  the  principal-stress  directions  coincide  with  the 

material-symmetry  direction.  It  must  be  remembered  that  the  principal  stresses 

and  their  directions  depend  directly  upon  the  original  stress  state 

(0,0,0)  with  respect  to  arbitrary  axes  x,  y . 
x y xy 

The  major  criticism  of  the  original  Ambartsumyan  model  insofar  as  it 
applies  to  filamentary  composite  materials,  is  that  it  does  not  relate  to  material- 
symmetry  directions,  which  are  of  major  importance  due  to  the  different  mechanisms 
of  tensile  vs.  compressive  stiffening  or  softening.  Thus,  the  model  does  not 
meet  criterion  1. 

The  Ambartsumyan  model  does  meet  criterion  2,  as  was  pointed  out  by 
Ambartsumyan  [9]  and  Tabaddor  [17]. 


Inspection  of  Eq.  (1) 
matrix  is  unsymmetrlc  when  0 


shows  that  the  Ambartsumyan-theory  compliance 
and  o 2 have  opposite  signs,  as  was  pointed  out 


matrix  is  unsymmetric  wnen  ui  ana  02  nave  opposite  signs,  as  was  pointea  ou 
by  Tabaddor  [l8]  and  Jones  [5].  Thus,  it  does  not  meet  criterion  3.  To  remedy 
this  disadvantage  of  the  Ambartsumyan  model,  the  variants  due  to  Isabekyan  and 
Khachatryan  [lOJ,  Shapiro  [ll],  and  Jones  [12]  were  introduced. 

‘4.  MODEL  II:  THE  ISABEKYAN-KHACHATRYAN  RESTRICTED-COMPLIANCE  MODEL 

This  model  was  introduced  by  Isabekyan  and  Khachatryan  [10]  and  is  a 
modification  of  the  Ambartsumyan  model  in  which  the  following  limits  are  imposed 
to  enforce  symmetry  of  the  compliance  matrix 

(-)  . v (+)  .V  (-) 


, <+> 

12 

, <+> 

361 
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, <+> 
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21 
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62 


Model  II  still  has  the  same  disadvantage  as  the  Ambartsumyan  theory  in 
that  it  does  not  meet  criterion  1.  In  fact,  the  limitations  imposed  upon  the 
coefficients  in  Eqs.  (3)  make  it  less  likely  to  meet  criterion  1.  Also,  it 
does  not  meet  criterion  2,  although  it  does  meet  criterion  3. 


5.  MODEL  III:  THE  SHAPIRO  FIRST-INVARIANT  MODEL 

This  model  was  introduced  by  Shapiro  [ll]  and  can  be  considered  to  be 
another  variation  of  the  Ambartsumyan  model,  Eqs.  (1),  in  which  now  the  bij  are 
determined  on  the  basis  of  the  sign  of  a combination  of  stresses  in  the  stress 
state,  namely,  the  first  invariant  of  stress  defined  as 


Oj  - (1/3)  (^  + o2  + a3) 


(1/3) (a  +0  + 0 ) / 
x y z 


For  the  plane-stress  case  considered  here, 

oT  - (l/3)(a1+a2)  = (l/3)(ax+a) 


(4) 


(5) 
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Thus,  here 


ij 


if  o > 0 


if  Oj  < 0 


(6) 


Again  the  major  disadvantage  of  Model  III  is  the  same  as  the  Ambartsumyan 
model  in  fail.ng  to  meet  criterion  1.  Also,  like  Model  II,  it  fails  to  meet 
criterion  2 but  does  meet  criterion  3. 


> •( 


6.  THE  JONES  WEIGHTED-COMPLIANCE  MODEL 

This  model  can  also  be  considered  to  be  a variant  of  the  Ambartsumyan 
model.  It  was  introduced  by  Jones  [12],  who  elaborated  upon  it  later  [5].  For 
the  plane-stress  case  [5], 


bu 


b21  = b12 


b .!♦> 

1J 

if 

°j  - 

b.^ 

IJ 

if 

°j  < 

b (+) 
b12 

if 

°1  > 

Vl2+> 

+ 

k2b12_) 

+ 

kb  (+) 
z°12 

if 

°1  < 

r 


11,  22,  61,  62 


if  o^  > 0 & o2  < 0 


(7) 


if  < 0 


o2  > 0 


where 


I 


i I + I °2 


I I + I °2 


(8) 


This  model  still  does  not  meet  criterion  1,  but  it  meets  both  criteria 


2 and  3. 


7.  MODEL  V:  THE  BERT  FIBER-GOVERNED  COMPLIANCE  MODEL 


This  model  was  introduced  by  Bert  [6]  and  differs  from  all  of  those 
previously  discussed  in  that  it  relates  stresses  and  strains  in  material-symmetry- 
axis  coordinates  (L,  T)  rather  than  principal-stress  coordinates.  Thus,  instead 
of  Eqs.  (1),  we  have  [6]: 


bLL  °L  + bLT  °T 
bLT  °L  + bTT  aT 
bss  °S 


(9) 
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where  subscript  S refers  to  shearing  action  relative  to  the  L,  T axes,  and 
r (+) 

bli  if  *fL  i ° 


b - < 
IJ 


^ 1£  °fL  ' 


°J 


f IJ  « LL,  LT,  TT 


(10) 


and  bgs  is  independent  of  stress.  Here  OfL  refers  to  the  fiber  stress  in 
the  fiber  direction.  A variety  of  approaches,  based  on  different  micromechanics 
assumptions,  can  be  used  to  approximate  the  critical  angle,  i.e.  the  angle  at 
which  an  applied  uniaxial  stress  system  must  be  oriented  (with  respect  to  the 
fiber  direction)  to  achieve  * 0 . Two  such  approaches  were  presented  in 

[6].  However,  it  is  not  necessary  to  know  this  angle  a priori;  the  test  data 
can  be  reduced  directly  in  such  a way  that  the  critical  angle  can  be  found 
directly  [6]. 


In  view  of  the  success  of  this  model  in  predicting  the  behavior  of  such 
a severely  bimodular  material  as  aramid-rubber  (see  Table  1) , it  is  believed 
that  this  model  comes  closer  to  meeting  criterion  1 than  any  of  the  other  models 
proposed  to  date.  Also,  it  meets  criteria  2 and  3. 

It  should  be  noted  that  the  Voigt  or  rule-of-mixtures  model  is  an 
excellent  micromechanics  model  for  representing  the  stresses  and  modulus  in  the 
fiber  direction  for  a single  layer  of  filamentary  composite  material.  This 
model  is  based  upon  the  assumption  that  the  fiber-direction  strain  in  the  fiber, 
matrix  and  composite  all  coincide.  This  affords  a simple,  practical  criterion 
(the  strain  in  the  fiber  direction)  to  determine  the  sign  of  the  fiber  stress 
at  any  given  point  in  a 'laminated-composite-material  structure. 


8.  FURTHER  DISCUSSION 

Three  other  unsymmetric  compliance  matrix  approaches  have  been  used  in 
reduction  of  experimental  data.  As  an  alternative  to  Model  V,  Bert  [6]  also 
applied  a nonsymmetric  monomodulus  model  in  reducing  aramid-rubber  data.  How- 
ever, he  demonstrated  that  this  model  was  not  as  satisfactory  as  Model  V.  Jones 
[19]  applied  both  the  original  Ambartsumyan  model  (Model  I)  and  an  unsymmetric 
orthotropic  bimodulus  modification  of  the  Jones-Nelson  nonlinear  model  [20] 
in  reducing  biaxially-stressed  ATJ-S  graphite  data.  He  also  specifically  studied 
the  significance  of  lack  of  compliance-matrix  symmetry  and  found  that  for  the 
material  considered  only  1 to  2%  difference  could  be  attributed  to  this  effect. 
Although  in  [5],  Jones  considered  compliance-matrix  symmetry  a necessary  criterion 
for  a consistent  material  model,  in  [19]  he  considered  the  unsymmetric  compliance 
matrix  to  be  "correct".  Regardless  of  the  validity  of  an  unsymmetric  compliance 
matrix,  since  it  makes  only  1 to  2 % difference,  it  is  clear  that  for  practical 
engineering  structural  computations,  a symmetric  matrix  is  to  be  desired. 


REFERENCES 

[ l]  A.E.H.  LOVE,  A Treatise  on  the  Mathematical  Theory  of  Elasticity. 
4th  Edn.  Dover,  New  York  (1944). 

[ 2]  L.R.G.  TRELOAR,  The  Physics  of  Rubber  Elasticity.  3rd  Edn.  Oxford, 
New  York  (1975) . 

[ 3]  F.P.  WOLF,  Polymer  13.  347  (1972). 


-9- 


[ 4]  B.  SAINT -VENANT , Notes  to  the  3rd  Edn.  of  Navier's  Resume  des  lemons  de 
la  resistance  des  corps  solides,  Paris,  175  (1864). 

[ 5]  R.M.  JONES,  AIAA  J.  15,  16  (1977). 

[ 6]  C.W.  BERT,  J.  Engng.  Matls.  Tech.,  Trans.  ASME  99H,  344  (1977). 

[ 7]  H.P.  PATEL,  J.L.  TURNER,  and  J.D.  WALTER,  Rubber  Chem.  Tech.  49,  1095 
(1976). 

[ 8]  S.A.  AMBARTSUMYAN,  Izv.  Akad.  Nauk  SSR,  Mekh.  no.  4,  77  (1965);  Engl, 
transl.  NTIS  Document  AD  675312  (1967). 

[ 9]  S.A.  AMBARTSUMYAN,  Mech.  Solids,  4 (3),  48  (1969). 

[10]  N.G.  ISABEKYAN  and  A. A.  KHACHATRYAN,  Izv.  Akad.  Nauk  Armyan.  SSR,  Mekh. 
22  (4),  25  (1966);  Engl,  transl.  NTIS  Document  AD  726600  (1971). 

[11]  G.S.  SHAPIRO,  Mech.  Solids.  1 (2),  85  (1966). 

[12]  R.M.  JONES,  AIAA  J.  9,  917  (1971). 

[13]  W.  VOIGT,  Lehrbuch  der  Kristallphysik,  p.  566.  Teubner,  Leipzig  (1928). 

[14]  R.A.  EUBANKS  and  E.  STERNBERG,  J.  Rat.  Mech.  Anal.  3,  89  (1954). 

[15]  B.M.  LEMPRIERE,  AIAA  J.  6,  2226  (1968). 

[16]  L.  BRUN,  J,  Elasticity  6,  89  (1976). 

[17]  F.  TABADDOR,  J.  Comp.  Matls.  3,  725  (1969). 

[18]  F.  TABADDOR,  AIAA  J.  10,  516  (1972). 

[19]  R.M.  JONES,  AIAA  J.  15.,  1436  (1977). 

[20]  R.M.  JONES  and  D.A.R.  NELSON,  JR.,  AIAA  J.  14,  709  (1976). 


{ 

I 


-10- 


MICROMECHAN ICS  OF  THE  DIFFERENT  ELASTIC  BEHAVIOR  OF 
FILAMENTARY  COMPOSITES  IN  TENSION  AND  COMPRESSION 


ABSTRACT 

It  has  been  known  for  a long  time  that  certain  kinds  of  fiber-reinforced 
composites  exhibit  quite  different  elastic  and  thermal-expansion  characteris- 
tics depending  upon  whether  the  fibers  are  In  tension  or  compression.  The  most 
dramatic  differences  have  been  noted  for  tire  cord/rubber  composites.  Two  dif- 
ferent mlcromechanistic  approaches  have  been  proposed  to  explain  these  observed 
differences  in  tension/compression  behavior:  the  mean-fiber-angle  approach  and 
the  elastically  supported  tie-bar/column  approach.  In  this  paper,  the  previous 
work  along  both  of  these  lines  Is  reviewed  and  extended. 


NOMENCLATURE 


- fiber  cross-sectional  area 

■ constant  of  proportionality  between  t and  d 

■ fiber  diameter 

* axial  Young's  modulus  of  the  fibers 

■ major  (fiber-direction)  Young's  modulus  of  the  composite 

■ axial  force  in  fiber 

* 8W,  Bq^o 

■ shear  modulus  of  the  matrix  material 

■ shear  modulus  of  the  composite  relative  to  1,  2 axes 

* area  moment  of  inertia  of  fiber  cross  section 

■ invariants  defined  In  equations  (37) 

■ foundation  modulus  of  matrix  material 

* fiber-direction  and  transverse  thermal  conductivities 

* half-wave  length  of  fiber  deflection  modal  shape 

■ bending  moment  acting  on  a fiber 

■ cos  8,  sin  8 

■ elastic  compliance  coefficients  (1,j  » 1,2,6) 

■ transverse  shear  force  acting  on  a fiber 

■ amplitudes  of  w and  w0,  respectively 

■ fiber  lateral  deflection  due  to  loading,  Initial  fiber  displace- 
ment from  a straight  line 

■ position  coordinates  parallel  and  normal  to  nominal  fiber  direction 

■ thermal -expansion  coefficient 

■ a of  fiber  and  matrix  constituents,  respectively 


This  part  of  the  report  Is  to  be  presented  at  the  Symposium  on  Mechanics  of 
Bimodulus  Materials,  sponsored  by  the  Applied  Mechanics  Olvlsion,  at  the 
ASME  Winter  Annual  Meeting,  New  York,  NY,  Dec.  1979,  and  will  be  published 
In  the  symposium  volume. 
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■ fiber  sinusoidal  parameter  and  Its  Initial  value 

■ temperature  change  from  the  strain-free  temperature 
• fiber-direction  and  transverse  strains,  total 

■ measured  fiber-direction  strain 

■ axial  normal  stress  on  fiber 

» Poisson's  ratio  of  the  matrix  material 

■ major  (axial-transverse)  and  transverse- transverse  Poisson's  ratios 
- angle  between  local  fiber  direction  and  the  nominal  fiber  direction 


Superscripts: 


indicates  that  this  Is  a property  of  the  wavy-flber  composite 
Indicates  that  this  Is  a mean  value 


Subscript: 

cr  denotes  critical  condition  at  which  onset  of  buckling  occurs 

INTRODUCTION 

In  1864,  Salnt-Venant  (l)1  presented  a mechanlcs-of  materials  analysis  of 
a beam  made  of  a material  having  a different  stress-strain  relation  when  load- 
ed In  tension  than  when  loaded  In  compression.  Although  this  phenomenon  has 
been  observed  in  certain  porous  cellular  media  such  as  sintered  metals  (2),  It 
has  been  observed  most  dramatically  In  fiber-reinforced  materials,  either  man- 
made (3.-6)  or  biological  (7.);  see  Table  1. 

A number  of  mechanistic  models  have  been  proposed  to  explain  the  different 
tension-compression  behavior  for  loading  In  the  fiber  direction.  All  of  these 
models  are  based  on  the  assumption  that  there  is  some  initial  curvature  in  the 
fibers,  so  that  they  are  straightened  and  thus  stiffened  by  tensile  loading 
and  curved  more  and  thus  made  less  stiff  by  compressive  loading.  All  of  the 
mechanistic  models  may  be  characterized  into  two  general  classes: 

1.  Tie-bar/column  on  elastic  foundation  model 

2.  Mean-fiber-angle  approach 


Table  1.  Some  Examples  of  Different  Elastic  Behavior 
in  Tension  And  Compression 


Investloator 

Ref. 

Material 

h!li 

Kotlarskl  & Karbasova 

(1) 

Fabric/rubber 

0.38 

Ducheyne  et  al. 

© 

Sintered,  porous 
stainless  steel 

0.1 

Zolotukhina  & Lepetov 

(i) 

Various 

fabric/rubber 

0.07  to  0.50 

Patel  el  al . 

(i> 

Polyester 
cord/rubber 
Aramld 
cord/ rubber 

0.017 

0.0034 

Ptfsfalvl 

(fi) 

Rayon 

cord/rubber 

0.0036 

Pearsall  & Roberts 

(D 

Myometrium 
(uterine  muscle) 

0.2 

In  the  ensuing  sections  both  of  these  models  are  described,  critiqued, 
and  extended. 


T 


Underlined  numbers  In  parentheses  designate  References  at  end  of  paper. 


METHOD  A:  TIE-BAR/COLUMN  ON  ELASTIC  FOUNDATION 
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Review  of  the  Literature 


Although  earlier  analyses  had  been  concerned  with  the  fiber  buckling  under 
compressive  loading  (8-TO),  perhaps  the  first  analysis  which  specifically  add- 
ressed the  different  fiber-direction  stiffnesses  in  tension  and  compression  was 
due  to  Herrmann  et  al.  (J1)  in  1967.  They  were  concerned  with  a relatively 
isolated  wire  in  a flexible  elastic  medium,  l.e.  a low-fiber-volume-fraction 
(Vf)  one.  The  application  with  which  they  were  concerned  was  a wire-reinforced 
solid-propellant  rocket  grain.  They  developed  two  different  analytical  pro- 
cedures, both  of  which  considered  the  matrix  as  a three-dimensional  elastic 
solid. 

Other  analyses  of  the  different  tension/compression  (T/C)  behavior  of 
filamentary  composites  have  been  concerned  with  more  closely  spaced  unidirec- 
tional composites.  Both  Bazant  (12)  and  Swift  (13)  used  two-dimensional  models 
analogous  to  those  used  by  Rosen  f£)  and  Chung  and  Testa  (14)  for  the  study  of 
fiber  buckling.  Of  course  there  are  two  distinct  fiber-buckling  modes:  ex- 
tensional  and  shear. 

It  has  been  shown  by  Lager  and  June  (15) , Greszczuk  (16),  Oavis  (17) , and 
Kulkarni  et  al . (18.)  among  others,  that  advanced  composites,  consisting  of 
boron,  graphite,  or  aramid  fibers  in  an  epoxy  matrix,  under  compressive  loading 
fall  in  the  shear  buckling  mode  rather  than  the  extensional  one.  However, 
these  same  composites  are  not  the  ones  which  exhibit  large  differences  in  T/C 
behavior  at  relatively  small  strains.  These  latter  ones  are  characterized  by 
low  Vf.  A simple  analysis  more  appropriate  to  the  latter  case  is  presented  In 
the  next  section. 

Before  leaving  the  topic  of  fiber  buckling  of  hlgh-Vf  composites.  It 
should  be  mentioned  that  Hayashi  (19),  Oavis  (J_7 ) > and  Wang  (20)  all  found  It 
necessary  to  take  into  account  the  effect  of  axial  compressive  stress  on  the 
shear  modulus  in  order  to  obtain  good  agreement  with  experimental  results.  . 
Also,  the  recent  work  of  Schaffers  (21_)  on  Intralaminar  fiber  buckling  under 
cylindrical  bending  should  be  mentioned. 

Recently  Baer  and  associates  (22-24)  have  conducted  extensive  research 
on  uncalclfied  collagen  tissues  and  have  developed  a mechanistic  model  based 
on  the  theory  of  an  extensible  elastica.  However,  their  model  does  not  take 
into  account  the  support  provided  by  the  surrounding  soft  tissue  (matrix 
material ). 


A Simple  Model 


The  analysis  presented  here  can  be  considered  as  an  extension  of  the  work 
of  (11 ) to  include  thermal -expansion/contraction  effects  due  to  temperature 
change,  with  a considerably  simplified  expression  for  the  foundation  modulus 
(8).  . 

The  equations  of  equilibrium  are: 


F'  • 0 
V a kw 

M'  - V + F (w0  + w)'  ■ 0 


Here  F?  axial  force  (i  for  tension/compression),  ka  elastic  foundation  modu- 
lus for  the  matrix.  M ? bending  moment,  V s transverse  shear  force,  ws  lateral 
deflection  of  fiber, w=  initial  deflection,  ( )'  s d(  )/dx,  and  x s axial 

position  coordinate. 

Assuming  the  fiber  to  be  a Bernoulll-Euler  linear  elastic  member,  the 
constitutive  relations  are: 

M - - EfI  w"  (4) 

F - AEf  [c  - («f  - oj  AT]  (5) 

Here  As  fiber  cross-sectional  area,  E^  s Young's  modulus  of  the  fiber  along 


' T / 


I 

1 


I 


t 


I 


Its  axis,  I * area  moment  of  Inertia  of  fiber  about  Its  centroid,  AT  * temper- 
ature change  measured  from  a strain-free  temperature,  Of  and  a fiber  and 
matrix  lineal  coefficients  of  thermal  expansion. 

In  view  of  equation  (1),  F Is  constant.  Combining  equations  (2),  (3), 
and  (4),  one  obtains  the  following  governing  equation: 

.IV 


Eflw 


F w"  ♦ kw  ■ F w0" 


(6) 


It  Is  very  likely  that  the  Initial  deflection  (w  ) Is  caused  by  a buckling 
of  the  fibers.  Thus,  It  Is  reasonable  to  assume  the  same  modal  shape  (dimen- 
sionless deflection  distribution)  for  the  Initial  deflection  and  deflection 
under  loading. 


W0  sin  e0x 


, w « W sin  bqx 


(7) 


Here  W and  W are  the  respective  amplitudes  of  w and  w , 8 a */i  , and  t * 
half-wave  length  of  deflection  modal  shape. 

Substituting  w and  w from  equations  (7)  Into  equation  (6),  one  obtains 
the  following  result: 

4 ♦ F B02  + k) 


W/WQ  " ' F B0V(EfV  + F Bo"  + k)  (8) 

Finally,  the  following  kinematic  relationship  Is  used  to  obtain  an  ex- 
pression relating  the  fiber-direction  strain  e,  to  the  deflection  amplitudes: 

t */fL 

el  " <8o/w)  | *'  dx 


Here 


f • (1/2)  C(w  *)2  - (w*)2]  + (F/AEf) 


no) 


tains: 


Substituting  equation  (10)  Into  equation  (9)  and  Integrating,  one  obtains 
c,  • (o/Ef)  + (am  - «f)  AT  + (B02/4)  (WQ2  - W2)  (11) 

Finally,  substituting  W/WQ  from  equation  (3)  Into  equation  (11),  one  ob- 


e1  * (o/Ef)  ♦ (<*m  * «f)AT 


+<W2)< 


(i-t— ^ 

l W* 


AEf[(o/Ef)  ♦ (om  ~ af)AT3Bfl 


AE^Cto/Ef)  + (om-<*f)AT]B0‘  + k 


(12) 


Special  Cases 

1.  No  Initial  curvature  (Bj,W  = 0): 
following  linear  relation: 


c1  « (o/Ef)  + (am  - Of)  AT 


Then  equation  (12)  reduces  to  the 


(12a) 


This  Is  the  linear  thermoelastlc  constitutive  relation.  Since  the  residual 
thermal  strain  (a  - af)  aT  Is  not  ususally  measured,  it  Is  customary  In  ten- 
sile tests  to  referencS  the  strain  from  the  value  of  zero  at  the  beginning  of 
the  test.  Thus,  the  measured  strain  Is 


c1(")  • o/Ef 


(12b) 


2.  High  tension  stress  (o  -v  large):  Then  c,  approaches  the  straight 
line  of  equation  (12a)  asymptotically  from  the  right. 

3.  Zero  stress  (osO):  Even  when  aT  a 0 also,  c.  Is  positive. 

4.  High  compression  stress  (-  c large):  The  limiting  case  Is 

Eflfl4  + F B2  ♦ k - 0 


mm 
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thus. 


- F • E*I$  ♦ k e 


-2 


03) 


Differentiating  equation  (13)  with  respect  to  a,  one  finds  the  Initial 
value  of  S to  be 


0cr  * (k/EfI)1/4 


(14) 


Since  I ■ wdy64  , where  d Is  the  fiber  diameter,  equation  (14)  Implies 
const. 


®cr  d 


or 


cr 


C d 


(15) 


This  result  Is  consistent  with  the  experiments  of  Rosen  (9)  and  the  elegant 
analysis  of  Sadowsky  et  al.  (25). 

To  evaluate  the  constant  C , we  apply  the  following  expression  for  the 
elastic  foundation  modulus  of  an  infinite,  isotropic  medium  surrounding  a rigid 
fiber  (8): 


16»Gm/[l  + 6(1  - 2vm)] 


For  a rubber  matrix,  v*  1/2. 

m 


Then 


k $ 16wG_ 


(16) 

(16a) 


Thus 


«cr  d - 2 ™ (W 


1/4 


(17) 


and 


<1/2 


I - «lcr  ■ «(Vfl 

5.  Decrease  in  temperature  only  (os  0,  aT  <0):  For  aT  ■ 0,  e,  > 0. 
Since  a > a,  for  practical  composites,  as  AT  Is  decreased  (aT  <0),  c,  also 
decreases.  The  limiting  case  Is  when 


(18) 


Efl04  + AEf  (am  - a^)62  AT  + k ■*  0 


Then 


|-AT|  - (EfI02  + k e‘2)/AEf  (offl  - «f) 


(19) 


This  has  a critical  value  of  e given  by  equation  (17)  and 

l-«lcr  ■ - ,f)  i 4(.ve f),/2/(^  - 

METHOD  B:  MEAN-FIBER-ANGLE  APPROACH 
Review  of  the  Literature 


(20) 


Although  an  earlier  analysis  was  concerned  with  Initially  wavy  layered 
media  (26).  the  first  analysis  specifically  addressed  to  the  use  of  this  ap- 
proach to  fiber- reinforced  composites  was  due  to  Tarnopol * skll  et  al.  (27). 
They  assumed  that  the  fibers  are  initially  sinusoidally  curved,  equation  (7a). 
Then  the  Inclination  at  an  arbitrary  axial  position  x is 


• • arctan  (e0WQ  cos  eQx) 


1 


(21) 
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Using  equation  (21)  In  the  well-known  transformation  equations  for  Young's  mod- 
ulus and  Integrating  over  a half-wave  length  t,  they  obtained  the  following 
result  for  the  fiber-direction  Initial  Young's  modulus  of  a wavy-flber-reln- 
forced  composite: 

(1/E,)  • (1/2E,)(2  ♦ f02)(l  ♦ f02)'3/2 

♦ (612_1  - 2v,2E1'1)(f02/2)(l  ♦ f02)’3/2 

♦ (1/2E2)[2  - (2  ♦ 3fQ2)(l  ♦ f02)"3/2]  (22) 

are  the  major  (fiber  direction)  and  minor  Young's  moduli  of  a 
composite  with  straight  fibers,  t.  Is  the  major  Young's  modulus 
of  the  wavy-fiber  composite,  f s B an<*  ^ vl2  are  the  resP*ct1ve 
major  shear  modulus  and  major  Poisson’s  ratio  *of  thestralght-flber  composite. 

For  f « 1,  the  following  approximate  expression  was  suggested  In  (27): 

Ej/E,  * [1  + (f02/2)(E1/G,2)]'1  (23) 

This  expression  Is  consistent  with  Bolotin's  result  (26). 

The  disadvantage  of  using  either  equation  (14)  or  (15)  Is  that  It  requires 
prior  knowledge  of  a straight-fiber  composite  having  the  same  fiber  volume 
fraction.  For  the  glass-fiber  reinforced  plastic  Investigated  experimentally 
In  (27).  It  was  found  empirically  that  an  Initial  composite  prestress  of  ap- 
proximately 10X  of  the  ultimate  tensile  strength  was  necessary  to  assure  that 
the  fibers  are  sufficiently  straight. 

In  an  alternate  approach  to  that  of  (27),  Nosarev  (28)  considered  the 
curved  fiber  as  a finite  number  of  straight^Tine  segements.  He  obtained  nu- 
merical results  which  predicted  reasonably  well  the  decrease  In  major  Young's 
modulus  and  increase  in  major  Poisson's  ratio  for  a copper-wire-reinforced 
epoxy  composite  as  the  fiber  curvature  Is  Increased. 

The  work  of  (26-28)  can  be  used  to  predict  onl/  the  Initial  modulus  (or 
Poisson's  ratio).  However,  In  certain  composites  with  a highly  flexible  ma- 
trix (such  as  rubber),  the  composite  modulus  may  change  quite  significantly 
with  loading  (especially  In  compression,  as  discussed  In  the  preceding  sec- 
tion). Thus,  the  analysis  of  Tabaddor  and  Chen  (29)  Is  quite  important  for 
such  composites,  since  it  permits  computation  of  the  complete  stress-strain 
curve  in  both  tension  and  compression  for  the  composite  In  the  fiber  direction. 
In  their  work,  they  replaced  f In  equation  (22)  with  f (s  6W)  and  developed 
a relationship  between  f and  f°. 

Makarov  and  Nikolaev  (30)  conducted  a series  of  experiments  on  aluminum- 
wire-  re  In  forced  natural  rubber  in  which  the  wires  were  provided  intentionally 
with  controlled  amounts  of  Initial  curvature.  The  ratio  of  Young's  moduli  of 
fiber  and  matrix  was  13,700  and  the  fiber  volume  fraction  was  0.06.  Jhey  found 
that  curvature  had  a negligible  effect  jjn  the  minor  Young's  modulus  (E-),  the 
transverse- transverse  Poisson's  ratio  (v--j,  and  the  transverse  thermaf-expan- 
slon  coefficient  (o2).  However,  for  their  test  conditions  they  measured  ap- 
proximately 502  decrease  in  E, , 48X  Increase  In  v,2,  and  15X  Increases  in  u,2 
and  a, . 

Vlshnevskll  and  Shlenskll  (31)  Investigated  the  effect  of  sinusoidally 
curved  fibers  on  transverse  thermal  conductivity  (Ic-)  of  the  composite.  They 
derived  an  expression  for  the  transverse  thermal  conductivity  which  Is  equiv- 
alent to  the  following  one  In  the  present  notation: 

*2  ■ (l/2)[k,f02  ♦ k2(2  - f02)]  (24) 

Here  k,  and  k.  are  the  conductivities  In  the  fiber  and  transverse  directions. 

It  should  be  mentioned  here  that  the  companion  equation  to  equation  (24) 
can  be  written  down  Inmedlately,  from  symmetry,  as  follows: 

♦ 0/2)[k1(2  - f02)  ♦ k2f02]  (25) 


Here  E,  and  E« 
unidirectional 


L 


i 


Tg'.  1 '■ ' 
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Also,  It  Is  noted  here  that.  In  view  of  the  fact  that  the  angular  transforma- 
tion equations  for  thermal  expansion  are  Identical  (second-rank  tensor)  to  those 
for  thermal  conductivity,  equations  (24)  and  (25)  are  also  applicable  to  ther- 
mal expansion  provided  that  the  k's  are  replaced  by  «'s. 

Tamopol'skll  et  al.  (32)  Investigated  the  effect  of  fiber  wavlness  In 
composites  with  anisotropic  fibers  (specifically  carbon  fibers).  They  found 
that  for  such  composites.  It  is  necessary  to  use  the  full  equation  (22)  rather 
than  the  simplified  equation  (23).  Van  Dreumel  and  Kamp  (33)  Investigated  ex- 
perimentally the  effects  of  fiber  waviness  on  both  tensile  strength  and  major 
Poisson's  ratio  of  carbon-fiber  reinforced  plastic. 

In  order  to  completely  characterize  a thin  sheet  material  for  use  In 
stress  analysis  and  structural  design,  one  needs  a complete  set  of  elastic 
stress-strain  relations,  for  fiber-direction  tension  and  compression  (34-35), 
not  just  the  Young's  moduli.  A rational  analytical  basis  for  such  a complete 
characterization  Is  developed  for  the  first  time  in  the  following  section. 

Complete  Characterization  of  Planar  Compliances 

A typical  repeating  element  Is  assumed  to  consist  of  a segment  of  a sin- 
usoidally curved  fiber  and  Its  surrounding  matrix  material.  The  fiber  curva- 
ture Is  assumed  to  be  planar  In  a single  plane  (xy)  with  a path  given  by 

y ■ Ma  sin  B.x  (26) 

Following  the  general  approach  used  by  Tabaddor  and  Chen  (29),  It  Is  as- 
sumed that  when  the  fiber  Is  subjected  to  a positive  (tensile)  strain  In  direc- 
tions 1 and  2,  respectively  parallel  and  normal  to  the  fiber  direction.  It 


takes  on  a new  sinusoidal  form  given  by 

y • W sin  sx  (27) 

where  the  relationship  between  s and  bq  and  between  W and  are  given  by 

(2v/B0)(l  + c,)  - 2*/fl  (28) 

W0(l  + c2)  - M (29) 

It  Is  assumed  that  shear  strain  has  no  effect  on  these  Inter-relationships. 
Combining  equations  (28)  and  (29),  one  obtains 

f » 6W  « B0W0  (1  + eg )/ ( 1 + «])  (3°) 

In  view  of  the  definitions  of  f and  of  Poisson's  ratio  (v.-2  " e^/e,), 

equation  (30)  can  be  rewritten  as  follows  for  the  uniaxial  case  only 

* f * fQ  (1  - c1)/(l  ♦ Cj)  (31) 


Equation  (31)  expresses  the  effect  of  c,  on  f(»  0W).  This  Is  a different 
approach  than  In  the  tie-bar/ column  one  (Method  A),  In  which  It  Is  assumed  that 
B » B.  and  W Is  found  by  solution  of  the  governing  differential  equation. 

"The  work  that  follows  is  based  on  the  elastic  compliance  rather  than  the 
planar  stiffnesses  or  the  engineering  moduli  (Young's  modulus,  Poisson's 
ratio).  There  are  several  reasons  for  this  choice: 

1.  Permits  direct  comparison  with  experimentally  measured  compliances. 

2.  Is  not  restricted  to  a particular  stress  state  (such  as  uniaxial,  bi- 
axial, or  even  triaxlal)  but  can  be  applied  directly  to  any  stress 
state. 

The  transformation  equations  relating  the  elastic  compliance  coefficients 
with  respect  to  axes  (r,  2),  oriented  at  an  angle  o with  respect  to  the  ortho- 
tropic material-symmetry  directions  (1,2),  to  those  with  respect  to  (1,  2) 
are  (36): 


U — M ■ ' 1 1 
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Here  m cos  t,  n i sin  e. 

In  writing  equations  (24),  the  contracted  Index  notation  was  used.  Sub- 
scripts 1 and  2 refer  to  normal  stress/straln  action  In  directions  1 and  2 re- 
spectively, while  subscript  6 refers  to  shear  stress/straln  action  related  to 
axes  1 and  2. 

To  average  the  properties  over  one  wavelength  of  fiber  waviness,  It  Is 
necessary  to  relate  direction  o to  the  curved-fiber  geometric  properties  W and 
e.  Working  toward  this  objective,  we  take  the  derivative  of  both  sides  of 
equation  (27): 


tan  e • y’  » WB  cos  ex 

Thus, 

m ■ cos  e ■ [l  + (we  cos  bx)2]"1^2 
n ■ sin  e ■ we  cos  ex/[i  ♦ (we  cos  ex)2]1^2 


(33) 

(34) 


are 


The  average  values  of  the  functions  of  m and  n appearing  In  equations  (32) 

mS  2 f L/4  m4  dx/2*  [ L/4  dx  » [1  ♦ (f2/2)](l  ♦ f2)  0/2 
J -l/4  / > -L/4 


(35) 


-L/4 

(TfV1'2 

m^fT  ■ mn3  » 0 

It  Is  now  possible  to  express  the  curved-fiber  composite  compliances 
ms  of  the  straight-fiber  unidirectional  compliances  as  follows:  J 
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It  Is  noted  that  since  m^n  > mn3  ■ 0,  the  wavy-flber  composite  remains 
orthotropic.  Use  of  equations  (30),  (35),  and  (36)  In  an  Iterative  manner  per- 
mits development  of  complete  nonlinear  strain  versus  stress  relationships  for 
the  composite. 
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oiscussion  and  concluding  remarks 


In  both  0/  the  mechanistic  analytic*!  methods  presented  here,  no  pro- 
vision was  made  for  transverse  shear  deformation  of  the  fibers.  In  fact,  there 
Is  some  controversy  as  to  which  Is  the  best  way  to  Incorporate  such  action;  see 
the  work  of  Roze  and  Klntsls  (37).  However,  a suitable  one  of  the  approaches 
described  In  (37)  could  be  Incorporated  Into  Method  A,  If  desired.  It  Is  noted 
that  In  Method  8 flexural,  muchover  shear,  properties  do  not  enter  explicitly 
Into  the  equations. 

The  experiments  of  Patel  et  al.  (5)  showed  that  the  difference  between 
tension  and  compression  properties  were~much  more  pronounced  for  rubber  re- 
inforced with  aramld  cord  than  when  reinforced  with  polyester  cord.  This  Is 
probably  related  to  the  very  low  compressive  strength  of  the  arainld  fibers 
themselves,  which  Greenwood  and  Rose  (38)  suggested  Is  due  to  separation  of  the 
microfibrils  comprising  the  fibers. 

Moncunlll  de  Ferran  and  Harris  (39)  found  experimentally  that  the  buckling 
mode  of  steel  wire/polyester  was  he! leaf  rather  than  planar  as  assumed  In  all 
known  analyses.  It  has  been  suggested  that  the  more  pronounced  difference  be- 
tween tensile  and  compressive  behavior  exhibited  In  the  tire-cord/rubber  com- 
posites of  (5-6)  as  compared  to  that  of  ordinary  wire-reinforced  rubber  (8) 
may  be  due  to  The  inherent  twisting  action  present  in  tire  cord,  which  ha?  a 
helical  geometry  (40).  However,  In  (5),  steel -cord/rubber  was  also  Investiga- 
ted and  no  appreciable  difference  between  the  tensile  and  compressive  stress- 
strain  behavior  was  noted. 

From  the  review  and  methods  presented  here,  It  can  be  concluded  that  the 
field  of  micromechanics  of  fiber-reinforced  composites  with  different  proper- 
ties In  tension  and  compression  Is  In  Its  Infancy.  Considerably  more  work  needs 
to  be  done,  among  which  the  following  appear  to  be  the  most  pressing: 

1.  Definitive  experiments  to  determine  which  of  the  two  methods  (A  and 
8}  presented  here  are  the  most  successful  In  predicting  composite 
material  behavior  (not  only  In  uniaxial  loading  but  also  under  multi- 
axial  conditions)  and  to  help  refine  or  even  combine  them. 

2.  Experiments  to  determine  the  role  of  fiber-matrix  Interfacial  adhe- 
sion on  the  gross  behavior. 

3.  Fracture  mechanics  analyses  and  carefully  controlled  experiments  with 
controlled  flaws. 

4.  Analyses  Incorporating  more  realistic  models  of  the  helical -strand 
action  In  tire  cords  (40). 

5.  Simple  means  of  Incorporating  the  different  T/C  behavior  Into  static 
and  dynamic  analyses  of  structural  components  such  as  plates  and 
shells  (41). 
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